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ABSTRACT 


The navigation and guidance analysis of a Mars 
probe launched from a manned flyby spacecraft is 
presented in this paper. This study was initiated with 
the probe and spacecraft separation at the Mars 
sphere of influence and was terminated with the probe 
arrival at a specified entry altitude and the spacecraft 
arrival at Mars periapsis. The disclosures of the 
study indicate that the results are decidedly dependent 
on the initial spacecraft errors at the time of separa- 
tion and on the system model errors assumed when 
processing navigation data. 



NAVIGATION AND GUIDANCE ANALYSIS OF A MARS PROBE 

LAUNCHED FROM A MANNED FLYBY SPACECRAFT 

By Thomas B. Murtagh, Flora B. Lowes, 
and Victor R. Bond 
Manned Spacecraft Center 

SUMMARY 


The navigation and guidance analysis of a Mars probe possibly to be launched 
from a manned spacecraft during Mars flyby missions in 1975 and 1976 is presented. 
The study is initiated with the probe and spacecraft separation at the Mars sphere of 
influence and is terminated with the probe arrival at a specified entry altitude and the 
spacecraft arrival at Mars periapsis. The relative state between the probe and space- 
craft is assumed to be perfectly known so that the uncertainty covariance matrix for the 
probe is equivalent to that of the spacecraft. The spacecraft position and velocity un- 
certainties are reduced by processing onboard sextant or theodolite measurements with 
a Kalman filter. The results of the study indicate that spacecraft periapsis -altitude 
dispersions between 1 and 6 n. mi. can be obtained with one or two midcourse velocity 
corrections totaling between 10 and 150 fps depending on (1) Mars radius uncertainty; 

(2) instrument accuracy; (3) sighting frequency; (4) type, number, and timing of cor- 
rections; and (5) the accuracy of their execution. The probe entry flight -path -angle 
dispersions ranged from 0. 13° to 0. 21° for the same range in velocity -correction- 
vector magnitude and for the same parameter variations as indicated for the spacecraft. 


INTRODUCTION 


Consideration is being given to the possibility of launching a Mars probe from a 
manned spacecraft during a Mars flyby mission between 1975 and 1980. Such a probe 
could be manned or unmanned. Basically, the probe must have the ability to land on 
the Martian surface and, after a stay time consistent with mission constraints, to 
launch and to rendezvous with the flyby spacecraft. 

The design of such a probe is influenced by many factors, not the least of which 
is the fuel required to complete the mission safely. The fuel required is, in turn, a 
function of many other parameters, such as (1) navigation system accuracy; (2) type 
and number of navigation measurements processed; (3) dispersions and uncertainties 
in the spacecraft trajectory at the time of separation; (4) type, number, and timing of 
midcourse velocity corrections and the accuracy of their execution; (5) distance from 
Mars when the probe is separated from the spacecraft; and (6) leadtime of the probe 
landing with respect to spacecraft periapsis passage. 



These last two parameters primarily affect the separation velocity of the probe 
with respect to the spacecraft. For a given leadtime, the separation velocity is 
approximately an inverse linear function of the distance from the planet; for a specified 
separation distance, the separation velocity is proportional to the leadtime required. 

To simplify the analysis of the probe-and-spacecraft navigation and guidance sys- 
tems, the separation distance and the leadtime were fixed so that the probe would sep- 
arate from the spacecraft at the Mars sphere of influence (SOI) and land on the Martian 
surface 1 hour ahead of the spacecraft periapsis passage. The relative state between 
the probe and the spacecraft was assumed to be known perfectly so that the uncertainty 
covariance matrices computed for the probe and the spacecraft were equal. The space- 
craft position and velocity uncertainties were reduced by onboard processing of sextant 
or theodolite optical measurements using a Kalman filter. The study was initiated with 
the probe and spacecraft separation and was terminated with the probe arrival at a 
specified entry altitude and the spacecraft arrival at Mars periapsis. 

The purpose of the study is to present the midcourse root-mean-square (rms) 
velocity corrections required for the probe to attain certain flight-path-angle disper- 
sions at entry (with corresponding vacuum periapsis-altitude dispersions) and for the 
spacecraft to attain reasonable altitude dispersions at the periapsis of the flyby hyper- 
bola. In the first section of the paper, the assumptions made in the study are 
discussed, followed by a section dealing with the reference trajectories used in the 
analysis. Next, a description of the navigation and guidance systems of the probe and 
flyby spacecraft is presented with a subsequent discussion of the digital computer pro- 
gram used to generate the results. The results and discussion section presents the 
navigation and guidance data in separate sections, further subdivided into discussions 
of either the probe or the spacecraft navigation and guidance systems. The pertinent 
mathematical developments are relegated to appropriate appendixes. 

The authors gratefully acknowledge the efforts of L. S. Cicolani of the NASA 
Ames Research Center in the generation of the covariance matrices presented in 
table I in this paper and wish to thank B. J. Garland of the NASA Manned Spacecraft 
Center for permission to reproduce figure 1 and for suggesting the launch dates and 
trip times used. 


SYMBOLS 


A(t) 

a 


B 

c(oX 2 ), s(aX 2 ) 


transformation matrix 
semimajor axis 

6 -by -3 matrix defined by equation (C5) 

Battin's auxiliary functions (The values of X are defined in eqs. (B25) 
and (B26). ) 


c 


ratio of planet-radius error to the radius of the observed planet 
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E, E(t) covariance matrix of state -vector uncertainty 

E E , Eq eccentric anomalies at entry and separation, respectively 

e eccentricity 

f true anomaly angle 

G(t) 6 -by- 6 guidance matrix 

Gj(t), GgCt) 3-by-3 submatrices of guidance matrix 

H, H (t) sensitivity matrix 

H„, Hq hyperbolic equivalent of eccentric anomaly at entry and separation, 

respectively 

h altitude 

I identity matrix 

i inclination angle 

K(t) weighting matrix 

L(t) 3 -by -3 matrix defined by equation (C8) 

!L (t) scalar defined by equation (C9) 

M(t) matrix defined by equation (All) 

N(t) covariance matrix of velocity-correction error 

P(t) matrix defined by equation (CIO) 

p semilatus rectum 

Q the measurable, as defined in equation (Al) (Q = /3 for sextant meas- 

urement) 

^2 2 2 

R = vx +y + z in equation (A6) 

/~2 2 

R' = vx + y in equation (A6) 

R(t) covariance matrix of measurement errors 
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position vector 
planet radius 

state vector (position and velocity) 

nominal time of arrival at target (or fixed terminal time) 
nominal time of periapsis passage 

time 

unit vector in the direction of ( ) 
down- range speed 

velocity vector 

matrix defined by equation (A8) 

vector defined by equation (D20) 

covariance matrix of state-vector dispersion 

components of the vector from the vehicle to the observed body 

gradient of periapsis radius with respect to position 

gradient of periapsis radius with respect to velocity 

gradient vector defined by equation (E6) 

gradient of flight -path angle with respect to position 

gradient of flight -path angle with respect to velocity 

reciprocal of semimajor axis 
right ascension 

star-planet included angle (sextant) 

celestial latitude and longitude of observed planet as viewed from 
spacecraft (theodolite) 

flight-path angle 


I 



Av 

Ay 

6 o 
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? 

v 2 
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velocity-correction vector 

probe entry flight-path-angle dispersion 

declination 

mean-squared engine-cutoff error 

mean- squared error in thrust-vector magnitude 

mean- squared error in thrust- vector orientation 
angle defined by equation (A12) 

gravitational constant multiplied by the mass of the central body 
parameter defined by equation (B14) 

sensor bias variance 
sensor variance 

total variance of observation errors 

state -transition matrix evaluated between t and t^ 

3-by-3 submatrices of the state-transition matrix 
right ascension of ascending node 


Subscripts: 

d 

E 

F 

h 


dispersion 

entry 

fixed time of arrival 

normal to the trajectory plane (in eq. (B9)) 


nom 


nominal 



p 

planet 

R 

required 

r o 

initial position 

s 

star 

s 

state (eq. (Al)) 

T 

terminal time 

u 

uncertainty 

V 

variable time of arrival 

vh 

vehicle-horizon 

vp 

vehicle-planet 

w 

trajectory plane 

6 

estimate 

7 T 

periapsis 

Superscripts: 

T 

transpose 

-1 

inverse 

+ 

value after measurement or correction 

mm 

value before measurement or correction 


Operators: 

6( ) variational operator (small deviation from reference value) 

£ ( ) expected value 

V ( ) gradient with respect to s 

s 
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ANALYSIS 


Assumptions 

The following ground rules are postulated for the analysis presented in this paper. 

1. The probe and the spacecraft are separated at the Mars SOI, with the probe 
landing approximately 1 hour prior to spacecraft passage through periapsis. 

2. Two sets of representative spacecraft uncertainty and dispersion covariance 
matrices at the Mars SOI were generated by propagating and updating an Earth-injection 
covariance matrix with rms position-and-velocity errors of 4 n. mi. and 16 fps. The 
first set, presented in table 1(a), used only Earth-based radar range and range-rate 
tracking data to update the spacecraft uncertainty matrix, and the second set, presented 
in table 1(b), used a combination of Earth-based and onboard navigation data to update 
the spacecraft uncertainty matrix during the Earth-to-Mars portion of the mission. A 
more detailed description of the systems used to generate these matrices can be found 
in reference 1. 

3. The reference trajectories used were assumed to be ideal conics, and the 
state-transition matrix used to propagate the errors was derived analytically for two- 
body conic trajectories. 

4. The relative state between the probe and the spacecraft is perfectly known so 
that the uncertainty covariance matrices of the probe and spacecraft are equivalent. 

The spacecraft covariance matrix is updated by processing onboard sextant or theodo- 
lite measurements with a Kalman filter. 

5. The total-variance equation of the observation errors is written as 


2 2 , . n ,2 2 

ct t = Uj + (c tan Q) + ct b 


(1) 


where cr^ is the standard deviation of the instrument error (equals 5 arc seconds) and 
is the standard deviation of instrument bias (equals 60 arc seconds). The param- 
eter c is defined as the ratio of planet-radius uncertainty to the planet radius; the val- 
ues used in this paper were 0. 001 and 0. 01 (corresponding to Mars-radius errors of 
2 and 20 n. mi. , respectively). The derivation and discussion of this equation are 
given in appendix A. 

6. Fixed time of arrival (FTA) guidance and variable time of arrival (VTA) 
guidance were used to compute the rms velocity corrections and appropriate target 
dispersions for both the probe and the spacecraft. No attempt was made to optimize 
the velocity-correction schedule for either of these types of guidance. 

7. A summary of navigation and guidance errors is presented in table II. 
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The Trajectory Characteristics 


The study is performed for the two Mars flyby launch opportunities of Septem- 
ber 20, 1975, and April 17, 1976. These trajectories will be designated as the short- 
long and the long- short cases, respectively. The trajectory characteristics are 
summarized in table in, and their heliocentric trajectories in the ecliptic plane are 
presented in figures 1(a) and 1(b). These trajectories have very nearly the same total 
trip times (671. 93 days and 680. 58 days, respectively). The trajectories are chosen 
primarily because they require near-minimum velocities at injection from the Earth 
(15 150 fps and 14 050 fps, respectively). The injection velocity is equal to the differ- 
ence between the perigee velocity of the escape hyperbola and the circular velocity at 
perigee. The perigee altitude at injection was 262 n. mi. for both cases. 

Both of these trajectories have hyperbolic asymptotes at Earth departure with 
declinations that are within +35°. It is pointed out in reference 2 that declinations in 
this range are necessary for good orbit determination with Deep Space Instrumentation 
Facility tracking stations. This range of declinations also permits departure trajec- 
tories with inclinations that are within the launch -azimuth (from Kennedy Space Center) 
safety constraints. 

It is important to consider the advantage of designing the flyby-probe mission so 
that it can be either a short-long or a long-short mission. The opportunities for 
minimum injection velocities for short-long (or long-short) missions are approximately 
765 days apart, but a long-short opportunity occurs 210 days after each short -long 
opportunity. Therefore, the long-short opportunity can be considered as a backup or 
alternate launch date to the short-long launch date. 

The initial conditions for the flyby spacecraft at the Mars SOI were obtained by 
matched conic techniques. These initial conditions are presented for the short-long 
and long-short cases in table IV. It is interesting to note the differences between these 
two trajectories within the Mars SOI. For the short-long case, the spacecraft is within 
the Mars SOI for approximately 38 hours, has a retrograde inclination with respect to 
the Martian equator of approximately 144°, and has a periapsis altitude of approxi- 
mately 106 n. mi. located in the northern hemisphere. For the long-short case, the 
spacecraft is within the Mars SOI for approximately 30 hours, has an inclination of 13°, 
and has a periapsis altitude of approximately 112 n. mi. located in the southern hemi- 
sphere. 

The probe trajectory is computed as shown in appendix B. The initial position 
vector is the same as that of the flyby spacecraft at the entrance to the Mars SOI. 
Specifying this initial position vector, along with the entry speed, the entry flight-path 
angle, the entry altitude, and the inclination, completely determines the probe trajec- 
tory, except for an ambiguity in the right ascension of the ascending node. The 
inclination is chosen to be the same as that of the flyby spacecraft to avoid plane change. 
The right ascension of the ascending node is also chosen to be the same as that of the 
flyby spacecraft to maintain the minimum distance between the flyby spacecraft and the 
probe. For this study, the probe entry altitude and flight-path angle assumed were 
315 000 feet and -20°. For the short-long and long-short trajectories, the entry speeds 
were 33 600 fps and 40 000 fps, respectively. 
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Description of Navigation and Guidance Systems 

The pertinent navigation- system equations are developed in appendix A along with 
the sensitivity vector, or matrix, which relates measurement deviations to state-vector 
deviations for the Mars-star included-angle measurement (sextant) and the spacecraft- 
Mars direction-of-line-of-sight measurement (theodolite). The geometry of the sextant 
included-angle measurement is illustrated in figure 2 ; the theodolite line-of-sight- 
direction measurement geometry is presented in reference^. The appropriate guidance- 
system error equations are developed in appendix C, and the derivation of the guidance 
matrix for FTA and VTA guidance is presented in appendix D. Appendix E is an outline 
of the procedure for computing the rms flight-path-angle dispersions or uncertainties 
and the radius-of-periapsis dispersions or uncertainties from the state dispersions or 
uncertainties. 


Description of Simulation 

A digital simulation program developed by the authors was used to generate the 
results presented in this paper. The basic components of this program are a control 
routine which generates trajectory and covariance-matrix time histories for both the 
probe and the spacecraft, a set of subroutines which compute the conic state-transition 
matrix for propagating errors and integrating the state vector along the conic, and 
another set of subroutines which update the covariance matrices as a result of a navi- 
gation measurement or a guidance maneuver. The rms position-and-velocity errors, 
computed from the square root of the trace of the covariance matrix, are presented in 
a locally level coordinate system which displays both in-plane and out-of-plane errors. 
The x-axis (altitude) of this coordinate system is along the radius vector to the probe 
or the spacecraft, the y-axis (range) is in the direction of the velocity, and the z-axis 
(track) is along the orbital angular-momentum vector. The errors in this system are 
designated as altitude, range, and track errors, and the time rates of change. 


DISCUSSION OF RESULTS 


Navigation Results 

Uncertainties in the estimate of the state of a vehicle as determined by navigation 
measurements represent a lower bound for the state dispersions. It is evident from 
equation (C4) in appendix C that state dispersions are never better than the uncertain- 
ties in the estimate. By investigating different navigation types to reduce the estimated 
uncertainties, the dispersions can also be reduced. Thus, as the navigation is im- 
proved, so is the ability to guide a spacecraft to a specified target. 

For navigation in the Mars-probe study, onboard optical determinations of an 
angle were used. Consideration was given to two optical instruments which might be 
used — the sextant and the theodolite. Although for parallel studies the theodolite 
tends to give somewhat better information than the sextant, the theodolite information 
was not found to be sufficiently better to justify presenting the theodolite data when the 
impracticability of the theodolite was considered. Thus, all data presented pertain to 
the star -planet included-angle measurement assumed to be made with an onboard 
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sextant and processed by a Kalman filter. Mars was the planet used for all measure- 
ments. However, the star used for each measurement was randomly chosen from a 
limited catalog of stars included in the simulation program. No attempt was made to 
optimize the choice of stars for the measurements. 

The effect of the previously described navigation method in reducing state -position 
uncertainties was investigated. This effect is presented for both of the trajectories 
described in table m, and in each case, the observation errors of table II and the initial 
uncertainty matrices of table I were used. Since the instrument bias a B (which is in- 
cluded in eq. (1)) can be related to the state, it is possible to incorporate in the 

state array and to have the filter act on it. Through this action, the bias could be 
reduced; consequently, the effect of the bias on the estimated uncertainties would be 
lessened. Only the instrument bias, not the basic instrument error, was estimated by 
this technique. The purpose of this estimation was to attempt to model errors which 
contribute to the overall error in estimating the state. 

Since the position uncertainty of principal interest is at periapsis for the space- 
craft and at vacuum periapsis for the probe, the results of the analysis are presented 
for those terminal points only. 

From previous studies, it has been found that a measurement frequency of less 
than 30-minute intervals does not significantly affect the periapsis position-error curve 
profile. Thus, the results presented are for measurement intervals of 30 minutes. 
Results of the navigation analysis are presented in figures 3 and 4. 

Spacecraft. - Figure 3 shows the spacecraft rms altitude uncertainties at Mars 
periapsis. These values are plotted against the time to Mars periapsis to illustrate 
approximately how well the periapsis-altitude uncertainty is known at any time along the 
trajectory. 

Figures 3(a) to 3(d) each contain four curves. The two solid-line curves repre- 
sent the rms periapsis-altitude uncertainty for a planet-radius error of approximately 
2 n. mi. (or c = 0. 001), and the two broken-line curves represent the rms periapsis- 
altitude uncertainty for a planet-radius error of approximately 20 n. mi. (or c = 0. 01). 
Of the two curves for each specified radius error, one curve represents the effect of the 
inclusion of instrument bias in the error model, and the other curve represents the 
effect of no instrument bias. The curves are marked accordingly. 

Figures 3(a) and 3(b) pertain to the short -long trajectory for two different sets of 
initial errors. The initial errors E^t Q ) and X^t^ , represented in figure 3(a), were 

generated from Earth-based radar tracking only. Those errors represented in fig- 
ure 3(b) were generated from the combination of Earth -based radar and onboard track- 
ing. It is evident from these two figures that the combination of Earth-based radar and 
onboard tracking gives lower initial uncertainties at the Mars SOI. Because of these 
lower initial uncertainties, those uncertainties computed for the time from the Mars 
SOI to Mars periapsis produce initially lower and continually smoother curve profiles 
than those produced by the initial errors generated from radar tracking alone. How- 
ever, the curves of both figures tend to converge to approximately the same four re- 
spective terminal values. 
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Figures 3(c) and 3(d) show the same curves as in figures 3(a) and 3(b) when gen- 
erated for the long-short trajectory. Practically the same conclusions can be drawn 
from these two figures as from figures 3(a) and 3(b). The slight differences in the pro- 
files of the curves and in the initial uncertainties are the result of the shorter time span 
of the trajectory between the Mars SOI and the Mars periapsis. 

For the respective case of interest, the effect of the calculation of spacecraft un- 
certainties can be determined from any of these four figures when a small planet-radius 
error of 2 n. mi. or a large planet-radius error of 20 n. mi. is considered. Also evi- 
dent in the figures is the effect of including an instrument bias in the error model with 
either of the two assumed radius errors. For the first few hours of navigation meas- 
urements, the instrument bias has the effect of keeping the uncertainty curve profiles 
from being lowered; this is followed by a general smoothing-off effect which is a result 
of the filter in the state estimation of the bias and of the tendency of the filter to de- 
crease the initial included bias value. It is possible to estimate the planet-radius er- 
ror in a similar manner, but this was not included in the study. Because of the relation 
between the estimated uncertainties and the state dispersions, which was previously 
pointed out, the study of the uncertainty curve profiles of figure 3 can aid in determining 
correction times. 

Probe . - Figure 4 pertains to the probe and presents the rms altitude uncertain- 
ties at vacuum periapsis for the short-long and long- short trajectories, respectively, 
and, also, the two types of initial errors in each case. The parts of figure 4 (a, b, c, 
and d) are presented in the same manner and sequence as the parts of figure 3 so that 
conclusions concerning the varying conditions may be drawn for the probe at vacuum 
periapsis in the same manner as for the spacecraft as the flyby-hyperbola periapsis. 

The initial projected probe-altitude uncertainties are slightly higher in value than 
the corresponding initial projected spacecraft uncertainties because of the degradation 
at separation in the covariance matrix E (tq)p R Q BE which results from an imperfect 

velocity correction execution. This degradation is evident from equation (C3). The 
degrading effect is noticeable in the probe uncertainty curves for the complete time 
span. However, the profiles of the curves of figures 3 and 4 are so similar that it is 
possible to draw approximately the same conclusions for the probe as for the space- 
craft. 


Guidance Results 

Spacecraft . - A representative set of spacecraft guidance data is summarized in 
table V for the short-long and long- short trajectories and for the two sets of initial un- 
certainty and dispersion matrices presented in table I. In the data presented, it was 
assumed that biased sextant-navigation measurements were processed every 30 min- 
utes with a 2-n. mi. Mars-radius error. Table V is included to illustrate the possible 
advantages of using VTA guidance as opposed to using FTA guidance and to compare the 
effect on total AV when making one or two midcourse velocity corrections for either 
of these guidance types. The target point for each type of spacecraft guidance is the 
periapsis of the flyby hyperbola. In FTA guidance, a velocity correction is computed 
such that target position-vector deviations are nulled; the VTA technique (referred to 
as radius-of -periapsis guidance in ref. 4) allows the down-range position error at the 
target to be free while minimizing the magnitude of the calculated correction. 
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The choice of the timing of the corrections presented was based upon the compar- 
ison of the rms estimate of AV computed from equation (C8) and the rms uncertainty 
in the estimate AV £ calculated from the square root of the trace of equation (CIO). 

This procedure, discussed in reference 5, is based on the assumption that a guidance 
maneuver will be commanded only after a certain amount of confidence in the trajectory 
estimate has been obtained. This assumption implies that AV^ should be smaller than 

AV, and the criterion chosen for this study was that AV g j£V < 0. 5. The philosophy 

of this choice of correction time will be discussed in more detail in the following sec- 
tion. 


The procedure for computing the radius-of-periapsis dispersion which results 
from a guidance maneuver is presented in appendix E. The significance of this particu- 
lar dispersion is that it is a measure of the "safe" planet periapsis passage of the 
spacecraft on the flyby hyperbola. 

To compare FTA and VTA guidance for the spacecraft, consider table V(a). As- 
sume that a single FTA correction is applied 0. 86 hour before periapsis, then the re- 
quired AV is 135 fps, with a resulting radius-of-periapsis dispersion of 1. 48 n. mi. ; 
a single VTA correction at the same time requires a AV of 117 fps, with a corre- 
sponding 1. 42-n. mi. radius-of-periapsis dispersion. The reduction in AV of 18 fps 
is obtained at the expense of an increase in down-range position error of 35 n. mi. This 
increase in down-range error is not significant, however, since this component of po- 
sition error can be mapped into a passage -of -periapsis timing error of a few seconds. 

Again, consider table V(a) to compare the effect on total AV of making either 
one or two velocity corrections to achieve a specified radius-of-periapsis dispersion. 
Assume that a dispersion of 1. 4 n. mi. is desired. A single FTA correction to produce 
this target dispersion requires a AV of 135 fps; if two corrections are applied, how- 
ever, the same dispersion results for a total AV of 17 fps. Using two corrections 
also enables the target dispersion to be further reduced to 1. 2 n. mi. for a total 
AV of 25 fps (FTA) and for a total AV of 23 fps (VTA). 

The preceding discussion involves data that were computed using initial uncer- 
tainty and dispersion matrices generated from only Earth -based radar tracking during 
the Earth-to-Mars portion of the mission. A comparison of tables V(a) and V(b) illus- 
trates the effect on the spacecraft midcourse velocity corrections in the Mars SOI of 
adding some onboard observations to the measurement schedule during the trans-Mars 
mission phase. The AV values in table V(b) are reduced, primarily because E^t^ 

generated from a combination of Earth -based radar and onboard tracking is much 
smaller than E^t^ generated from only Earth-based tracking. The initial dispersion 

matrix X(t Q ) is approximately the same for both types of tracking. 

The preceding discussion for tables V(a) and V(b) (short -long trajectory) applies 
equally well to tables V(c) and V(d) (long-short trajectory). 

Probe . - Figure 5 presents a comparison of the rms estimate of AV, the uncer- 
tainty in the estimate AV € , and the probe entry flight- path-angle dispersion Ay, when 
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biased sextant-navigation measurements are processed at 30-minute intervals with a 
2-n. mi. Mars-radius uncertainty. Figures 5(a) and 5(b) are plots for the short-long 
trajectory and the two sets of initial covariance matrices, while figures 5(c) and 5(d) 
are similar plots for the long-short trajectory. These plots illustrate the philosophy 
which might be used to select the time of a single velocity correction to achieve a spec- 
ified target Ay dispersion while satisfying the criterion that AV £ y'AV < 0. 5. 

Consider figure 5(a) and assume that an entry flight -path -angle dispersion of 0. 16° 
is desired. A single correction 5 hours prior to the arrival at the nominal entry alti- 
tude will produce a Ay dispersion of 0. 16 ° for a AV of 20 fps. The uncertainty in 
the velocity estimate at this point is 2 fps, so that AV^/aV = 0. 1 and the selection 

criterion is satisfied. Now consider figure 5(b) and assume a required Ay dispersion 
at entry of 0. 155°. This dispersion can be achieved as early as 5 hours prior to entry, 
and the Ay plot is flat from that time on. Therefore, to keep the AV as small as 
possible, the earliest correction time is selected (i. e. , 5 hours before entry). The 
rms AV is 3 fps, but AV £ is 2 fps, so that AV^AV = 0. 67 and the selection cri- 
terion is not satisfied. Similar discussion would apply to figures 5(c) and 5(d). 

A representative set of probe guidance data is summarized in table VI. For 
these data, biased sextant- navigation measurements were processed at 30-minute in- 
tervals, assuming a*2-n. mi. Mars-radius uncertainty. The target of the FTA guid- 
ance technique is the probe vacuum periapsis, and this technique attempts to null the 
position-vector deviations at that point. The target of the VTA guidance technique is 
the nominal entry altitude, which allows the down -range position error to be free 
while nulling the entry flight -path -angle deviation. The procedures for calculating the 
probe radius-of-vacuum-periapsis dispersion and the entry flight -path -angle dispersion 
are outlined in appendix E. The Ay dispersions presented in table VI are computed 
for fixed radius (e. g. , at the nominal entry altitude). The initial uncertainty and dis- 
persion matrices for the probe were computed from the E^q) and X^t ^ of the space- 
craft at the time of separation and were properly degraded as a result of the separa- 
tion maneuver. 

The data in table VI indicate that there is very little apparent difference between 
the FTA and VTA guidance for the probe. This might be expected to be the case if 
equation (D16) is compared to equation (D31). Equation (D16) presents the submatrices 
G^(t) and Gg(t) of the guidance matrix G(t) for FTA guidance; equation (D31) repre- 
sents these submatrices for VTA guidance nulling Ay at the target. The submatrix 
Gg(t) is the same for both guidance schemes, and G^(t) for the VTA guidance reduces 

to Gj (t) for FTA guidance if 


0 


12 1 - T - 2 T ( < ^ ) 21 ~ ^ 22^12 
(Zi T + — 2 T< ^22^ > 12 _1 )— T 



<<C ^12 ^11 
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The interpretation of the data presented in table VI follows the same pattern as 
the interpretation of the spacecraft guidance data in table V. For example, consider 
table Vl(a). A single FTA correction executed 0. 83 hour before entry requires a AV 
of 111 fps with a corresponding Ay of 0. 18°; if two corrections are applied, with the 
second correction applied 0. 83 hour before entry, the total AV is 42 fps with a Ay 
of 0. 16°. Again, the use of two corrections has produced a significant saving in AV. 

If VTA guidance is used for these same two corrections, the total AV is 44 fps with a 
Ay of 0.15°. The difference to be noted is that while the total AV is approximately 
the same, the Ay error is better controlled with VTA guidance, whereas the radius- 
of-vacuum-periapsis error is better controlled with FTA guidance. This is to be 
expected because of the nature of the two guidance laws being compared. The radius- 
of-vacuum-periapsis dispersion is a measure of the entry corridor, and by multiplying 
that number by some appropriate "safety" factor, the corridor width attainable by the 
guidance system can be computed. Normally, this safety factor is 6, so that for the cases 
just cited, the corridor widths would be 8. 64 n. mi. and 9. 18 n. mi. , respectively. 

Tables VT(b) and VI(d) present only single- correction results because two correc- 
tions to the probe trajectory were not necessary in these cases. These data resulted 
from the smaller initial probe uncertainty covariance matrix used. The discussion of 
the data in table VI(c) would be similar to that presented for table VI(a). 


CONCLUDING REMARKS 


The navigation and guidance analysis of a Mars probe launched from a manned 
flyby spacecraft has been presented. The study was initiated with the probe and space- 
craft separation at the Mars sphere of influence and was terminated with the probe ar- 
rival at a specified entry altitude and with the spacecraft arrival at Mars periapsis. 

The results of the .study indicate that spacecraft radius-of-periapsis dispersions 
between 1 and 6 n. mi. can be obtained for a total velocity between 10 and 150 fps. The 
probe entry flight -path -angle dispersions varied from 0. 13° to 0. 21° (corresponding to 
an entry corridor of 10 n. mi. ) for approximately the same range in velocity. These 
results are very dependent on the initial spacecraft errors at the time of separation and 
on the system model errors assumed when processing navigation data. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, January 10, 1968 
981-30-89-00-72 
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TABLE I. - SPACECRAFT UNCERTAINTY AND DISPERSION COVARIANCE MATRICES AT MARS SPHERE OF 
INFLUENCE IN MARS EQUATORIAL COORDINATE SYSTEM, ft, fps 
(a) Generated from Earth- based radar tracking only 



+0. 10795743E+11 

+0. 65967692E+09 

+0. 33814298E+11 

+0. 18949483E+04 

-0. 37702034E+03 

+0. 10354719E+04 



+0. 18766645E+10 

+0. 70484404E+10 

-0. 53478838E+03 

+0. 10425965E+04 

+0. 42295007E+03 




+0. 11982285E+12 

+0. 39853662E+04 

+0. 16151761E+04 

+0. 43253749E+04 

II 




+0. 17066236E-02 

-0. 45197054E-03 

+0. 26474908E-03 



(Symmetric) 



+0. 16697713E-02 

+0. 20050646E-03 


— 





+0. 21117166E-02 


+0. 73819030E+12 

+0. 15401839E+12 

-0. 93253504E+11 

-0. 86159 568E+07 

-0. 18965936E+07 

+0. 15584321E+07 



+0. 17476617E+12 

-0. 70977560E+11 

-0. 18317749E+07 

-0. 20811097E+07 

+0. 93663495E+06 




+0. 28377198E+12 

+0. 15426505E+07 

+0. 95201726E+06 

-0. 18607330E+07 

o 

ti 




+0. 10291022E+03 

+0. 22766192E+02 

-0. 18675098E+02 



(Symmetric) 



+0. 25288863E+02 

-0. 11509438E+02 







+0. 22327910E+02 


cji 



TABLE I.- SPACECRAFT UNCERTAINTY AND DISPERSION COVARIANCE MATRICES AT MARS SPHERE OF 


INFLUENCE IN MARS EQUATORIAL COORDINATE SYSTEM, ft, fps - Concluded 
(b) Generated from Earth- based radar and onboard tracking 


+0. 59 298541E+09 

-0. 49410256E+09 
+0. 64494082E+09 

+0.34391141E+09 
+0. 16991 21 7E+09 
+0. 15493640E+10 


(Symmetric) 


"+0. 73238785E+12 

+0.15312751E+12 
+0. 17436720E+12 

-0. 11271050E+12 
-0.74533639E+11 
+0. 21538386E+12 


(Symmetric) 



+0. 33825529E+03 

-0. 31423879E+03 

-0. 47339011E+02 

-0. 28329428E+03 

+0. 43679336E+03 

+0. 61234597E+02 

-0. 17214476E+01 

+0. 65811396E+02 

+0. 15218244E+03 

+0. 16175887E+00 

+0. 25352728E-01 

-0. 18663014E-01 


+0. 37405348E-01 

-0. 11504193E-01 



+0. 62743453E-01 

-0. 11453944E+08 

-0. 23871738E+07 

+0. 17817307E+07" 

-0. 24086431E+07 

-0. 26909 167E+07 

+0. 11632760E+07 

+0. 17784900E+07 

+0. 115759 28E+07 

-0. 31164389E+07 

+0. 179448 20E+03 

+0. 37444451E+02 

-0. 28045894E+02 


+0. 41732681E+02 

-0. 18035724E+02 


+0. 51504544E+02 



TABLE II. - NOMINAL ROOT -ME AN- SQUARE ERROR VALUES, la 


Errors in making velocity correction: 

Magnitude, percent 

Direction, deg 

Cutoff, fps 

Observation noise: 

Instrument error, a^, arc sec . . 
Radius uncertainty/planet radius, c 
Instrument bias, a_, arc sec . . , 

Jo 


1 

1 

0. 5 


5 

0 . 01 / 0.001 

60 


TABLE IB. - COMPARISON OF SHORT-LONG (SEPTEMBER 20, 1975, LAUNCH) 
AND LONG-SHORT (APRIL 17, 1976, LAUNCH) MARS FLYBY 
TRAJECTORY CHARACTERISTICS 


Trajectory designation 

Short- long 

Long- short 

Julian date of launch from Earth 

2 442 675.0 

2 442 885.0 

Earth injection- velocity 

magnitude, fps 

15 150 

14 050 

Declination of Earth departure 

asymptote, deg 

33.31 

-30.24 

Outbound-trip time, days 

133.29 

582.21 

Return-trip time, days 

538.64 

98.37 

Time in Mars SOI, hr 

37.74 

29.92 

Periapsis altitude at Mars, n. mi 

106.09 

112.67 

General location of periapsis at Mars . . . 

Northern hemisphere 

Southern hemisphere 

Inclination to Mars equator, deg 

144.21 

13. 10 

Lighting conditions 

Periapsis -dark 

Periapsis -light 

Entry velocity at Earth, fps 

47 900 

51 300 
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TABLE IV. - INITIAL POSITION AND VELOCITY OF FLYBY SPACECRAFT AT MARS SPHERE 


OF INFLUENCE IN MARS EQUATORIAL COORDINATE SYSTEM FOR SHORT- LONG 
AND LONG- SHORT TRAJECTORIES, n. mi., fps 


Short- long 

X 

-0. 30318650E+06 

Y 

-0. 57007078E+05 

Z 

+0. 52589475E+05 

X 

+0. 27000551E+05 

i 

Y | 

+0. 52456193E+04 

• 

Z 

-0. 45743130E+04 

Long- short 

! 

X 

+0. 30693083E+06 

Y 

+0. 5883609 2E+05 

Z 

-0. 16437793E+05 

X 

! -0. 34589177E+05 

• 

Y 

1 -0. 639043 58E+04 

• 

Z 

+0. 17966413E+04 


.! 



TABLE V. - SPACECRAFT ROOT-MEAN-SQUARE AV SUMMARY WITH BIASED SEXTANT-NAVIGATION 
MEASUREMENTS PROCESSED EVERY 30 MINUTES, ASSUMING A 2-NAUTICAL MILE MARS-RADIUS ERROR 
(a) Short-long trajectory; E(t Q ), X generated using Earth-based radar tracking only 


Guidance 

scheme 

! 

Root- 

mean- 

square 

AVj, 

fps 

Time of AV^ 

measured from 
Mars periapsis, 
hr 

Root- 

mean- 

square 

av 2 , 

fps 

Time of AV 2 

measured from 
Mars periapsis, 
hr 

Total AV, 
fps 

Radius- of- 
periapsis 
dispersion, 
n. mi. 


10. 46 

11.36 

- 

i 

10. 46 

5.38 

Fixed 

15.15 

7.86 

- 

- 

15.15 

3.42 


41.39 

2.86 

- 

- 

41.39 

1.86 

time 

135.00 

.86 

- 

- 

135. 00 

1.48 

of 

10. 46 

11.36 

1.62 

6.86 

12.08 

3.41 

arrival 

10. 46 

11.36 

5.37 

2.36 

15.83 

1.80 


10. 46 

11.36 

6.85 

1.86 

17.31 

1.49 


10. 46 

11.36 

14.80 

.86 

25.26 

1.25 


9. 31 

11.36 

- 

. 

9.31 

5.44 

Variable 

13. 47 

7. 86 

- 

- 

13.47 

3.45 


36. 57 

2. 86 

- 

- 

36.57 

1.83 

time 

117.65 

.86 

- 

- 

117.65 

1.42 

of 

9. 31 

11.36 

1.55 

6.86 

10. 86 

3.48 

arrival 

9. 31 

11.36 

5.16 

2.36 

14. 47 

1.83 

9.31 

11.36 

6. 58 

1.86 

15. 89 

1.51 


9.31 

11.36 

14.08 

.86 

23. 39 

1.26 
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TABLE V. - SPACECRAFT ROOT-MEAN-SQUARE AV SUMMARY WITH BIASED SEXTANT -NAVIGATION 
MEASUREMENTS PROCESSED EVERY 30 MINUTES, ASSUMING A 2- NAUTICAL MILE MARS- RADIUS ERROR - Continued 


(b) Short-long trajectory; 


X M 


generated using Earth-based radar and onboard tracking 


Guidance 

scheme 

I 

Root- 

mean- 

square 

AV p 

fps 

Time of AVj 

measured from 
Mars periapsis, 
hr 

Root- 

mean- 

square 

AV 2 , 

fps 

Time of AV 2 

measured from 
Mars periapsis, 
hr 

— 

Total AV, 
fps 

Radius- of - 
periapsis 
dispersion, 
n. mi. 


5.82 

11.36 

- 

- 

5. 82 

3.92 

Fixed 

8.46 

7. 86 

- 

- 

8.46 

2.83 


23.49 

2.86 

- 

- 

23.49 

1.63 

time 

79.25 

.86 

i 

- 

79.25 

1.28 

of 

5. 82 

11.36 

1.05 

6.86 

6. 87 

3.00 

arrival 

5. 82 

11.36 

3.69 

2.36 

9.51 

1.68 


' 5.82 

i 11.36 

4.74 

1,86 

10. 56 

1.42 


5.82 

11.36 

10.28 

.86 

16.10 

1.21 


! 3.51 

11.36 

j 

- 

3.51 

4.23 

Variable 

5.12 

, 7.86 

- 

- 

5.12 

3.02 


13.98 

! 2.86 

- 

- 

13.98 

1.63 

time 

45.05 

.86 

i 

- 

45. 05 

1.21 

of 

3.51 

i 11.36 

; 1.06 

6.86 

4. 57 

3.28 

arrival 

' 3.51 

i 11.36 

1 

| 3.86 

2.36 

7.37 

1.84 


3.51 

j 

| 11.36 

4. 96 

1.86 

i 

8.47 

1.53 


' 3.51 

11.36 

10. 70 

.86 

14.21 

1.28 




TABLE V.- SPACECRAFT ROOT- MEAN-SQUARE AV SUMMARY WITH BIASED SEXTANT-NAVIGATION 
MEASUREMENTS PROCESSED EVERY 30 MINUTES, ASSUMING A 2-NAUTICAL MILE MARS- RADIUS ERROR - Continued 
(c) Long- short trajectory; E(t^, X(t Q ) generated using Earth-based radar tracking only 


Guidance 

scheme 

Root- 

mean- 

square 

AV r 

fps 

Time of AV^ 

measured from 
Mars periapsis, 
hr 

Root- 

mean- 

square 

av 2 , 

fps 

Time of AV 2 

measured from 
Mars periapsis, 
hr 

Total AV, 
fps 

Radius- of - 
periapsis 
dispersion, 
n. mi. 


19.73 

7.46 

- 

- 

19.73 

3.88 

Fixed 

27.00 

5. 46 

- 

- 

27.00 

2. 59 

time 

49.85 

2.96 

- 

- 

49.85 

2.11 

of 

153.84 

.96 

- 

- 

153.84 

1. 58 

arrival 

19.73 

7.46 

1.73 

4.96 

21.46 

2.61 

19.73 

7.46 

3.37 

2.96 

23.10 

2.12 


19.73 

7. 46 

5. 30 

1.96 

25. 03 

1.51 


19.73 

7.46 

22.86 

.46 

42.59 

1.21 


15.18 

7.46 

- 

- 

15.18 

3.87 

Variable 

20.75 

5.46 

- 

- 

20.75 

2. 53 

time 

38.13 

2.96 

- 

- 

38.13 

2. 01 

of 

115.76 

.96 


- 

115.76 

1.41 

arrival 

15. 18 

7.46 

1.60 

4.96 

16.78 

2.61 

15.18 

7.46 

3. 18 

2.96 

18.36 

2.12 


15.18 

7.46 

4.98 

1.96 

20.16 

1.50 


15. 18 

7.46 

21.63 

.46 

36.81 

1.20 


to 
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TABLE V. - SPACECRAFT ROOT- MEAN- SQUARE AV SUMMARY WITH BIASED SEXTANT- NAVIGATION 


MEASUREMENTS PROCESSED EVERY 30 MINUTES, ASSUMING A 2-NAUTICAL MILE MARS- RADIUS ERROR - Concluded 


(d) Long- short trajectory; E^t^, X ( t o) 


generated using Earth-based radar and onboard tracking 


Guidance 

scheme 

Root- 

mean- 

square 

* v i> 

fps 

Time of AV^ 

measured from 
Mars periapsis, 
hr 

Root- 
mean- 
s qua re 
AV 2 , 

fps 

Time of AV 2 

measured from 
Mars periapsis, 
hr 

Total AV, 
fps 

Radius -of- | 

periapsis 
dispersion, 
n, mi. 

| 

14.83 

7.46 

- 

- 

14. 83 

3.13 

Fixed 

20.32 

5. 46 

“ 

- 

20.32 

2.27 j 

i 

* time 

37.65 

2.96 

- 


37.65 

1.87 

1 

i of 

117.58 

.96 

- 


117.58 

1.38 j 

arrival 

14.83 

7.46 

1.18 

4.96 

16.01 

2.40 ’ 


14.83 

7.46 

2.40 

2.96 

17.23 

1.97 | 


14.83 

7.46 

3.84 

; 1.96 

18.67 

1.45 

i 

14.83 

7. 46 

16.83 

1 * 46 

31.66 

1.17 


8.40 

7.46 

- 

- 

8.40 

3.28 

Variable 

11. 50 

5. 46 

- 

- 

11. 50 

2.34 

time 

21.16 

2.96 

- 

- 

21.16 

1.81 

of 

64.39 

.96 

- 

- 

1 64.39 

1.20 

arrival 

8.40 

7. 46 

1.11 

4.96 

9.51 

i 

2. 58 1 


8.40 

7.46 

i 2.34 

2.96 

10.74 

2.13 


8.40 

7.46 

' 3.75 

1.96 

12.15 

1. 55 


8.40 

: 7.46 

1 16.63 

.46 

| 25.03 

1.23 





TABLE VI.- PROBE ROOT-MEAN-SQUARE AV SUMMARY WITH BIASED SEXTANT-NAVIGATION 
MEASUREMENTS PROCESSED EVERY 30 MINUTES, ASSUMING A 2-NAUTICAL MILE 

MARS -RADIUS ERROR 

(a) Short-long trajectory; E^t^, X^t Q j generated using Earth -based radar and onboard tracking. 


Guidance 

scheme 

Root- 

mean- 

square 

AV r 

fps 

, Time of AV^ 

measured from 
entry, 
hr 

Root- 

mean- 

square 

AV 2 , 

fps 

— 

Time of AV 2 

measured from 
entry, 
hr 

Total AV, 
fps 

Radius- of - 
vacuum- 
periapsis 
dispersion, 
n. mi. 

Entry 

flight-path- 

angle 

dispersion, 

deg 


7. 18 

13.33 

- 

- 

7. 18 

It. 41 

0.21 


10. 14 

9.83 

- 

- 

10. 14 

5.45 

.17 

Fixed 

25.79 

3.83 

_ 


25.79 

2.36 

.16 

time 

111.22 

.83 

- 

- 

111.22 

1.56 

. 18 

of 

7.18 

13.33 

3.34 

8.83 

10.52 

5.19 

.17 

arrival 

7.18 

13.33 

8.05 

3.83 

15.23 

2.38 

.15 


00 

t-H 

13.33 

35.29 

.83 

42.47 

1.44 

.16 


7. 18 

13.33 

81.24 

.33 

88.42 

1.40 

. 17 


7.20 

13.33 

- 

- 

7.20 

11.41 

0.21 


10. 18 

9.83 

- 

- 

10. 18 

5.45 

.17 

Variable 

26.05 

3.83 

_ 

_ 

26.05 

2.39 

. 15 

time 

116.47 

.83 

- 

- 

116.47 

2. 19 

. 15 

of 

7.20 

13.33 

3.35 

8.83 

10.55 

5.20 

.17 

arrival 

7.20 

13.33 

8.14 

3.83 

15.34 

2.39 

. 15 


7.20 

13.33 

37. 18 

.83 

44.38 

1.53 

.15 


7.20 

13.33 

93.36 

.33 

100. 56 

1.82 

. 15 


co 

co 




TABLE VI. - PROBE ROOT -ME AN- SQUARE AV SUMMARY WITH BIASED SEXTANT- NAVIGATION 
MEASUREMENTS PROCESSED EVERY 30 MINUTES, ASSUMING A 2-NAUTICAL MILE 

MARS- RADIUS ERROR - Continued 

(b) Short- long trajectory; E(tp), X(t^ generated using Earth-based radar and onboard tracking 


Guidance 

scheme 

Root- mean- 
square AV, 
fps 

Time of AV 
measured from 
entry, 
hr 

Radius- of- 
vacuum-periapsis 
dispersion, 
n. mi. 

Entry 

flight-path- angle 
dispersion, 
deg 

Fixed 

3.83 

3.83 

2.18 

0.15 

time 

5.20 

2.83 

1.97 

. 15 

of 

8.01 

1.83 

1.69 

. 15 

arrival 

40.07 

.33 

1.24 

| .16 

t 

Variable 

3.88 

3.83 

2. 19 

0.15 

time 

5.29 

2.83 

1.98 

. 15 

of 

8.24 

1.83 

1.70 

. 15 

i arrival 

1 

i 47.08 

.33 

1.37 

.15 


TABLE VI. - PROBE ROOT-MEAN-SQUARE AV SUMMARY WITH BIASED SEXTANT-NAVIGATION 
MEASUREMENTS PROCESSED EVERY 30 MINUTES, ASSUMING A 2-NAUTICAL MILE 

MARS- RADIUS ERROR - Continued 

(c) Long- short trajectory; E(t Q ), X(t Q ) generated using Earth-based radar tracking only 


Guidance 

scheme 

Root- 

mean- 

square 

AV r 

fps 

Time of AV^ 

measured from 
entry, 
hr 

Root- 

mean- 

square 

av 2 , 

fps 

\ 

Time of AVg 

measured from 
entry, 
hr 

Total AV, 
fps 

Radius- of- 
vacuum - 
periapsis 
dispersion, 
n. mi. 

Entry 

flight-path- 

angle 

dispersion, 

deg 


9.77 

9.87 


- 

9. 77 

7.25 

0. 19 


18.58 

5.37 

- 

- 

18.58 

3.47 

.14 

Fixed 

70.75 

1.37 

_ 


70. 75 

1.67 

.13 

time 

238. 92 

.37 

- 

- 

238. 92 

1.43 

.14 

of 

9. 77 

9.87 

5. 13 

5.37 

14. 90 

3.56 

. 14 

arrival 

9. 77 

9.87 

11.91 

2.37 

21.68 

2.06 

.13 


9.77 

9.87 

20.36 

1.37 

30. 13 

1.45 

.13 


9.77 

9.87 

69. 12 

.37 

78. 89 

1.32 

. 13 


9.80 

9.87 

- 

- 

9. 80 

7.25 

0. 19 

Variable 

18.68 

5.37 

- 

- 

18.68 

3.48 

. 14 

time 

72.22 

1.37 

- 

- 

72. 22 

1.72 

.13 

of 

256. 53 

.37 

- 

- 

256. 53 

1.92 

.13 

arrival 

9.80 

9.87 

5. 15 

5.37 

14.95 

3.56 

. 14 

9.80 

9.87 

12.02 

2.37 

21.82 

2.07 

.13 


9.80 

9.87 

20.70 

1.37 

30. 50 

1.47 

.13 


9.80 

9.87 

73.27 

.37 

83.07 

1.46 

.13 


to 

ui 



TABLE VI. - PROBE ROOT -ME AN- SQUARE AV SUMMARY WITH BIASED SEXTANT- NAVIGATION 
MEASUREMENTS PROCESSED EVERY 30 MINUTES, ASSUMING A 2-NAUTICAL MILE 

MARS- RADIUS ERROR - Concluded 

(d) Long- short trajectory; E generated using Earth-based radar and onboard tracking 


Guidance 

scheme 

Root-mean- 
square AV, 
fps 

Time of AV 
measured from 
entry, 
hr 

Radius- of - 
vacuum-periapsis 
dispersion, 
n. mi. 

Entry 

flight- path- angle 
dispersion, 
deg 

Fixed 

6.01 

2.37 

1.87 

0.13 

time 

10.39 

1.37 

1.47 

. 13 

of 

16.15 

.87 

1.34 

. 13 I 

I 

arrival 

35.51 

.37 

1.22 

. 13 

Variable 

6.01 

2.37 

1.88 

0.13 

time 

10. 52 

1.37 

1.48 

.13 

of 

16.46 

.87 

1.36 

. 13 

arrival 

37.23 

.37 

1.29 

.12 




(a) Launch date — September 20, 1975 (short-long). 

Figure 1. - Projection of free-return trajectories into the ecliptic plane. 


27 





28 


(b) Launch date — April 17, 1976 (long -short). 
Figure 1. - Concluded. 





rms altitude uncertainty at periapsis, n.mi. 



Time to Mars periapsis, hr 


(a) Short-long trajectory; E(t Q ), X(t Q ) generated from 
Earth-based radar tracking only. 

Figure 3. - Spacecraft rms periapsis -altitude uncertainties with sextant -navigation 
b measurements processed every 30 minutes. 


30 


Time to Mars periapsis, hr 


(b) Short -long trajectory; E ft^ , X^t^ generated from 
Earth-based radar and onboard tracking. 


Figure 3. - Continued. 
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Figure 3. - Continued. 
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(d) Long-short trajectory; E^t^, -^(^o) £ enera ted from Earth -based 

radar and onboard tracking. 


Figure 3. - Concluded. 
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(a) Short-long trajectory; E^t^, ^‘C'o) S enerated from 

Earth- based radar tracking only. 

Figure 4. - Probe rms vacuum periapsis -altitude uncertainties with 
sextant -navigation measurements processed every 30 minutes. 
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(b) Short -long trajectory; E (t^ , X^t^ generated from 
Earth-based radar and onboard tracking. 


Figure 4. - Continued. 
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(c) Long- short trajectory; E(t Q ), X(t Q ) generated from Earth-based 

radar tracking only. 
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Figure 4. - Continued. 
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(d) Long- short trajectory; E(t Q ), X (t Q ) generated from Earth-based 

radar and onboard tracking. 


Figure 4. - Concluded. 
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(a) Short-long trajectory; E^ 0 j, X^t Q ) generated from Earth- based radar tracking only. 

Figure 5. - Comparison of the rms estimate of velocity correction, the rms uncertainty in the 
estimate, and the resulting probe entry flight -path -angle dispersion. 
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(b) Short-long trajectory; generated from Earth- based radar and onboard tracking. 
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Figure 5. - Continued. 
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Figure 5. - Concluded. 


APPENDIX A 


NAVIGATION SYSTEM EQUATIONS 


The material developed in this appendix can generally be found in references 3 
to 5. For onboard spacecraft navigation during a Mars flyby mission, the principal 
types of measurements that can be made are optical determinations of an angle. In 
this study, two optical instruments — the sextant and the theodolite — were consid- 
ered. 


In simulating the use of the sextant, measurements were made of the star-planet 
included angle. For each observation, a star was chosen randomly from an abbrevi- 
ated catalog of stars incorporated in the simulation program. No attempt was made to 
optimize the selection of the star used. Since the study originates at and terminates 
in the Mars SOI, the observed planet was always Mars. Similarly, for the theodolite 
measurement, the planet Mars was always the observed body. 

For each measurement, it is necessary to compute the H matrix, that is, the 
matrix of partial derivatives which relate measurement deviations to state deviations. 
The H matrix is associated with each type of observation and is expressed by 


H(t) = V g Q 


(Al) 


where Q is the measurable. 


Sextant 

For the sextant measurement, Q is the star-planet included angle /?, and H(t) 
is a seven-dimensional vector expressed by 


H(t) T 




1 


(A2) 


where the 1 component is representative of the derivative of the observable with re- 
spect to an associated bias, the component = 0, and the component 3/3/Sr^ 

is derived by treating the star-planet included -angle measurement as a special case 
of the Sun-planet measurement (ref. 5). The geometry of the star-planet included- 
angle measurement is shown in figure 2. 


42 



Thus, assuming a vehicle- star vector r and a vehicle -planet vector 

— s 

the derivation, the equation for the included angle /3 can be expressed by 


r 

-vp 


for 


cos j3 = 



(A3) 


By taking first-order deviations of equation (A3), expanding the obtained form, 

and passing to the limit of the increasing r , the following equation for the star-planet 

s 

included- angle measurement can be obtained. 


6/3 = 





- /u • u \u 
\~s -vp) -vp 

| r vp| Sl ” P 



<A4) 


where u and u are the unit vectors from the spacecraft to the star and to the 
— s — vp 

planet, respectively. 

From equation (A4), 3(3/ dr is obtained; consequently, H(t) for the star- 

/ vp 

planet included-angle measurement is completed. 


Theodolite 


For the theodolite measurement, Q would be the measurables, celestial latitude 
and longitude, of the observed body (denoted by /3^ and ^ respectively), and H(t) 

would be an 8-by-2 matrix expressed by 


H(t) T 


dP 1 


3r 

dv 

-vp 

-vp 


v 2 

ar 

av 

-vp 

-vp 


(A5) 
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In this equation as in equation (A2), the derivatives of the observables with respect to 
the velocity and with respect to an associated bias would equal 0 and 1, respectively. 
The derivative with respect to the position is expressed by 


a/3j 


xz 

yz 

Sr 

-vp 


R 2 R’ 

R 2 R’ 



y 

X 

Sr 

L _vp J 


R' 2 



2 r,2 

z - R 


R 2 R' 



(A6) 


where x, y, and z are components of the vector from the vehicle to the observed 
body, R = Vx 2 + y 2 + z 2 , and R’ = Vx 2 + y 2 (appendix C of ref. 3). 

As a result of the navigation measurement and, consequently, of the evaluation 
of the sensitivity matrix H, knowledge of the vehicle state can be improved by using 
the information obtained to update the covariance matrix of state uncertainties E 
which is propagated along the trajectory by the equation 


E(t)=0(t,t o )E(t o )«(t,t o ) T 


(A7) 


where <p ^t, t^^ is the state -transition matrix evaluated between tg and t. To insure 

that E(t) remains at least positive semidefinite, E is replaced in the propagation 
procedure by a matrix W so that 


E(t) = W(t)W(t) T (A8) 


The replacement of E by W is used only in the propagation procedure and not in the 
updating procedure. 

When a measurement is made, the updated value of the E matrix is calculated 

from 


E(t) + = E(t)“ - K(t)H(t)E(t)“ 


(A9) 
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The weighting matrix K(t) is chosen so that the mean- squared position and velocity 
uncertainties are simultaneously minimized, 


K(t) = E(t)"H(t) T M(t)“ 1 


(A10) 


where 


M(t) = H(t)E(t)"H(t) T + R(t) 


(All) 


The covariance matrix of measurement errors R(t) is considered, for this study, to 
have four components: (1) error in the instrument, (2) error in the knowledge of the 

radius of the observed planet, (3) error associated with instrument bias, and (4) error 
caused by the uncertainty in the observed position of the planet. However, since the 
observed planet is the central body for the conic section, the last error component is 
zero (ref. 4), and consideration is given to only the first three error components. 

For the sextant measurements, R(t) is a scalar, and for the theodolite measure- 
ments, R(t) is a 2-by-2 matrix. 

The first component to be considered is the variance of the optical instrument 

2 

tfj , which is already known. The second error component, the lack of knowledge of 

the radius of the observed planet, is inversely proportional to the distance of the 
planet from the vehicle. Thus, the change in 9 (one -half the angle subtended by the 
observed planet) must be considered. 

Referring to figure 2 for the geometry used in the derivation of the second error 
component of R(t) 


r 

sin 9 ~ 2— (A12) 

r 

vp 


By taking deltas of the above equation and making the appropriate substitutions, 
the following equation is obtained. 


60 = tan 9 -^-2 (A13) 

P 
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from which the parameter c is defined as the ratio of the planet-radius error to the 
radius of the observed planet, or 


c 



(A14) 


The third component, error associated with instrument bias, is given an initial 
value and is then estimated in the filter along with the state variable. 

Thus, the total- variance equation of observation errors can be expressed as 


R(t) = a T 2 



2 , 2 a 
+ c tan 9 


+ CT r 


(A15) 
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APPENDIX B 


SOLUTION FOR THE CONIC TRAJECTORY THAT SATISFIES THE 
TERMINAL ENTRY CONDITIONS OF v , y , h E 


The material presented in this appendix is based in part on reference 5. It is 
assumed that the probe is at the initial position vector r^ and that it is desired that 

the terminal or entry speed v„, flight-path angle y„, and altitude h„, be satisfied. 

Hi Hi Hi 

The inclination i of the trajectory is specified. The geometry of this problem 
is shown in figure B-l. The velocity required v^v^ = Xprobe) at *0 to satis -ty 

these conditions will now be derived. 

The semimajor axis of the trajectory is computed from 



where r„ = h + r . If a > 0, the trajectory is elliptic; if a < 0, the trajectory is 

Hi ill P 

hyperbolic. 

The parameter p is found from 


P = 



2 . 2 
sin y E 


P 


(B2) 


The eccentricity of the trajectory is then 


e = V 1 " a (B3) 

The radius of periapsis is found from 

r^ = a(l - e) (B4) 


47 


I 



It should be noted that if the trajectory is elliptical, the constraint 


r. < 2a - r 
0 7T 


(B5) 


must be imposed. 

From the polar equation of the orbit, the angle 0 between r_ and r_ may be 
found from the true anomalies 



and 

by setting 

6 = f E ' f 0 

The normal to the trajectory plane is computed from 

u, = sin fi sin i u - cos sin i u + cos i u 
— h —x — y — z 


(B6) 


(B7) 


(B8) 


(B9) 


The inclination i is specified, and R is the right ascension of the ascending node 
which may be computed from r^ and i, except for an ambiguity which may be freely 

chosen. 

From figure B -2, the right ascension a q and declination 5^ of r^ maybe 
computed from 


a n 



(BIO) 
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and 



(Bll) 


Also from figure B-2 


ft = a 


0 


a 


(B12) 


for one choice of ft, or 


ft = <Xq + a + 7T 


(B13) 


for the other choice of ft. For both cases, a is computed from 


sin a = 


tan 6 q 
tan i 


(B14) 


Equation (B14) implies the limits on the inclination which must be specified to make 
the equation meaningful; that is 


tan i > tan 6 q 


(B15) 


From figure B-l, it is evident that the entry position and velocity may be 
written 


-E 


o -0 

= r^ cos 9 — + r„ sin 6 u 


‘ E 


E 


~w/ 


(B16) 
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and 


-E ~ V E 


-0 


cos ( e + >e) + V E sln ( e + 


(B17) 


where the unit vector lies in the trajectory plane and is defined by 


H w 'H h x !l I 


(B18) 


The velocity required may now be found simply from the triple crossproduct 

£o x (^o * — r) = — o x (— e * — e) (b19) 

which is true since r A x v_ = r„ x v„ is the angular- momentum vector which must 
remain constant. After considerable algebra, equation (B19) may be solved for v R 
as follows. 


In r v 

^R = F^ v 0 COSy 0 + ^-¥^ siny l 


The magnitude of v R is found from 



The initial flight -path angle y ft may be found from 


— 0 X — R 


-E 


X — E 


(B20) 


(B21) 


(B22) 
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or 


r 0 V R Sin y 0 = r E V E Sin y E 


y 0 


= sin 


■ l/ r E V E SlU y E \ 
\ r 0 V R I 


(B23) 


The time to entry may be computed from the equation 


t E = r ° ‘ X 2 c(ffX 2 ) + (1 - r 0 «)X 3 s(«X 2 ) + TqX (B24) 

yfi* 

where a = 1/a and c(aX 2 ) and s(«X 2 ) are Battin's auxiliary functions. The 
value X is found from 

( e e- e o) 

yfa' 

if a > 0, and from 

„ ( h e- h o) 

yj-a 

if a<0. If the trajectory is elliptical (a > 0), the quantities Eg and E Q are the 
eccentric anomalies and are found from 


(B25) 


(B26) 


tan 




- e 
e 



(B27) 
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If the trajectory is hyperbolic (a < 0), the quantities H E and Hq are found from 


tanhg) 



(B28) 


When these equations are applied to determine the v_ of the probe, the inclina- 

tion and right ascension of the ascending node are the same, or very nearly the same, 
as those of the spacecraft. 
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£ ^ is a unit vector out of page 


Figure B-l. - Geometry of trajectory which satisfies the entry 
conditions v £ , y E> r E = r p + hg. 



Figure B-2. - Geometry of the two trajectories possible that 
satisfy the entry conditions. 




APPENDIX C 


GUIDANCE SYSTEM ERROR EQUATIONS 


The material presented in this appendix may be found in references 4 and 5. To 
perform a complete statistical evaluation of a guided flight, it is necessary to propa- 
gate and update the covariance matrix of the state uncertainties E(t) and the covari- 
ance matrix of the state dispersions X(t). Both matrices are propagated using the 
equations 


E(t) = 0(t,t o )E(t o )0(t,t o ) T (Cl) 

and 

X(t) = 0 (t, t 0 ) X (t Q ) 0 (t, t Q ) T (C2) 


where 4>(t, t^j is the state -transition matrix evaluated between t^ and t. If a naviga- 
tion measurement is made at time t, the matrix X(t) remains unchanged. However, 
if a guidance correction is commanded, then both E(t) and X(t) are updated accord- 
ing to the following equations (ref. 5). 


E(t) + = E(t)~ +BN(t)B T 


(C3) 


and 


X(t) + = [ I + G(t)J [X(t)' - E(t)'] [ I + G(t)] T + E(t) + 
The 6 -by -3 matrix B is defined by 


B = 


(C4) 


(C5) 
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where I is a 3 -by -3 identity matrix. The 6 -by -6 guidance matrix G(t) is written 


G(t) = 


0 

G x (t) 


0 

G 2 (t) 


(C6) 


The derivations of the submatrices G^(t) and G 2 (t) for the analysis are presented in 
appendix D. 

The matrix N(t) is defined as the covariance matrix of the velocity -correction 
error and is derived and discussed in references 4 and 5. The difference between the 
two derivations is that in reference 5 it is assumed that the inaccuracy in establishing 
a commanded velocity vector is caused only by random errors in thrust-vector orien- 
tation and magnitude. Reference 4 includes a third random error, designated as an 
engine -cutoff error. In this study, N(t) is used as derived in reference 4 and may be 
written as 


N(t) = ? 2 L(t) 


+ 



[i (t)I " L(t)] + 


e 2 Ut) 


(C7) 


where 


L(t) = [GjW G 2 (t) 


[ X(t) - E(t)] [G l( t) G 2 (t]] T 


and 


£ (t) = trace L(t) 


(C8) 


(C9) 


2 2 2 

The parameters C , r? , and e are the mean-squared values of the errors in thrust 
magnitude, orientation, and cutoff, respectively. The rms estimate of the velocity 
correction is computed from equation (C8). The rms uncertainty in this estimate is 
computed from 


P(t) = [GjW 


G 2 (t)j E(t) [GjW G 2 (t)] T 


(CIO) 
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APPENDIX D 


DISCUSSION OF GUIDANCE FOR FIXED TIME OF ARRIVAL 
AND FOR VARIABLE TIME OF ARRIVAL 


A small velocity correction is to be made at time t to a trajectory to guarantee 
that certain terminal conditions will be met. Two general classes of terminal con- 
straints which may be considered are fixed time of arrival (FTA) constraints and var- 
iable time of arrival (VTA) constraints. This appendix is an interpretation of FTA and 
VTA guidance presented in references 6 and 7, but the final results are obtained in 
terms of the familiar state-transition matrix. A particular application of VTA guid- 
ance, the guidance-to- entry nulling flight- path- angle errors which are not discussed in 
these references, is also presented. 


At the instant before the velocity correction, the terminal position r^, is re- 
lated to the position and velocity just prior to correction time by 


-t" = -T ( D1 ) 


Lack of a subscript on r and v means that they are evaluated at correction time t; 
a subscript T means that they are evaluated at the terminal time T. 

Taking linear deviations of equation (Dl) gives 



6r 



9v 


(D2) 


At the instant after the velocity correction, the terminal position £ T + is related to the 
position and velocity immediately after correction time by 



CD3) 
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Taking linear deviations of equation (D3) gives 





(D4) 


Since there can be no instantaneous change in position, r + = r 


r, and 


6r + = 6r" = 6r 


(D5) 


If the correction is assumed to be small and the trajectories are linearly related 


and 


3£ t " 3r T + 


9r~ 

- + 
dr 

Sr-p 

0rrp + 



3v" 

3v + 


(D6) 


(D7) 


These partials are the submatrices 0^ and 0^ °f the state-transition matrix. 

(For the conic case, these are given in part in ref. 5. ) Equations (D2) and (D4) may 
be written, using these simplifications, as 


6r T " = 0 n 6r + 0 12 5 v" (D8) 

and 

6r T + = 0 1;l 6r + 0 12 6 v + (D9) 
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Note that 6v 


and 6v + are related to the nominal trajectory by 


6v 



-nom 


(DIO) 


and 


6v + = v + - v 
- — -nom 


(Dll) 


The velocity-correction vector Av is found from 


Av = v 


+ 


v 


\ 


or 




(D12a) 


Av = 6v + - Sv~ J 


Now, solving equation (D8) for 6v” and equation (D9) for 6v + , and by differencing the 
results, Av is obtained as 


Av = 




(D12b) 


Replacing 6r T by equation (D8), Ay becomes 


A - = ^12 _1 ( 6 -T + " ^n 6 £ " ^12 6 ^") (D13) 

FTA Guidance 

In this type of guidance, it is desired to null the position deviations 5r T + at the 
fixed terminal time T, that is 
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(D14) 



= 0 


Then Av becomes 


A -F = “^12 _1< ^11 6 - " 


or 




(D15) 


Ay F = G^Sr + GgSv 


where 


G 




-I 



(D16) 


VTA Guidance 

In this type of guidance, the terminal time T is allowed to be free. For ex- 
ample, let the component of the position deviation parallel to the velocity at T be 
given by 


5 £ t + = ~-T 6T (D17) 

Substituting into equation (D13), Ay for VTA is given by 

A ^v= ^la' 1 (-V T - ^n 5 I - *!2 5 i) O 518 ' 

The VTA correction can then be expressed in terms of the FTA correction as 

Ay v = Ay F - 0 12 -1 y T 6T (D19) 
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By defining 




equation (D19) becomes 


AVy= AVp, - W oT 


Since there is an additional degree of freedom remaining, 6T may be chosen, if de- 
sired, to minimize Av y . From equation (D21) 

Av y 2 = AVp, 2 - 2W- Av p 6T + W- W 6T 2 1 


= -2W • Av p + 2W • W 5T = 0 


W- Av_ 
— -F 

W-W 


Substituting back into equation (D21) 


/ W W\ 

A -V = l 1 " W • W / A -F = G l 6 - + G 2 6 - 
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where 


G 


1 


G 2 




W W T \ 

W-W j (p l2 ^11 


T 

W W 1 
W-W ' 


/ 


(D24) 


VTA Guidance, Null 6y T 

The extra degree of freedom in VTA guidance may also be used to control an ad- 
ditional terminal deviation if desired. As an example, the errors in flight -path angle 
at the terminal time may be nulled. 

The deviation in flight-path angle is shown in appendix E to be a function of posi- 
tion and velocity deviations at T such that 


^T=Zi T S£ T ++ Z 2 T 6v T + 


(D25) 


But 


Sr T + = 0 1;l 6r + 4> 12 6 v + 


(D9) 


and 


6 I T +=0 21^ + ^22^ + 


(D26) 


Solving these two equations for 6y T + and eliminating 6v + 


5 V = 0 21 6 I + ^ 22 ^ 12 ”* ( 6 V ' * ll 5 l) 


(D27a) 
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or 


5 X T += 21 ‘ <P 22 <P 12~ lct> ll) 8 * + ^12" l6 XT + ( D27b ) 


Substituting this equation into equation (D25) to eliminate 

6> T = (-1 + —2 ^22^12 ) 6 ~T + + ~2 ! {^2\ " ^22^12 ^ll) 6 ! (° 28 ) 

Using equation (D17), setting = 0, and solving for 6T yields 


V (*21- Vl2' Si)* 

6T = — (D29) 


This equation is the time deviation required for the deviation in flight -path angle 
to be nulled. Substituting equation (D29) into equation (D18) 


where 


(Av v ) = G. 6r + G 9 6 v" 

V ~ VJ 6y T = 0 1 “ 1 ~ 


(D30) 


G l = - 


^12 X - r-2 (^21 " ^ 22^12 ^ 11 ) 

(z i T + Z 2 T0 22 0 12 *)-T 


^‘^ll 


G 2 = - X 


(D31) 


These guidance laws are applied directly to the guidance-to-entry problem discussed 
in the body of this paper. The terminal conditions subscripted by T are the same as 
the entry conditions subscripted by E. 
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APPENDIX E 


COMPUTATION OF ROOT -MEAN -SQUARE FLIGHT -PATH -ANGLE 
AND PERIAPSIS -RADIUS ERRORS 


The basic equation which relates the flight -path angle y to the position r and 
the velocity v is 


r • v 


cos y - 


rv 


(El) 


from which 


6 y = Z 1 T 6r + Z^Sv 


(E2) 


If the computation of the flight -path -angle error is desired for fixed time, then 


r x (r x v) 


Zl = 2 

rlrxvl 


(E3) 


If the calculation requires flight -path -angle error for fixed radius, then 


— 1 ~ |r xy| 


(E4) 


In either case (fixed time or fixed radius), written 


—2 = 2 


(r x v) x v 


v r x v 


(E5) 
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T / T T\ 

Defining = I Z^ j, equation (E2) becomes 


where 


Therefore 


6y = 


Z T 5s 



2 T T 
8y = Z 5s6s Z 

Taking expected values of equation (E8) gives 

£ (dy 2 ) = Z T f(6s6s T )z 


(E6) 


(E7) 


(E8) 


(E9) 


Since X(t) is defined by ^fos5s T l the mean -squared flight -path -angle dispersion is 
computed from ' ' 


°yd = Z(t) T X(t)Z(t) 


(E10) 


Similarly, the mean-squared flight -path -angle uncertainty can be computed from 


V 2 = Z(t) T E(t)Z(t) 


(Ell) 


The derivation of the periapsis -radius errors can be computed by at least two methods. 
The first method is quite similar to that presented for the calculation of flight -path - 
angle errors and can be found in reference 4. The radius of periapsis can be 
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written as a function of the position and velocity at any point along the conic, from 
which 


6r w = yj 5r + y 2 6v (E12) 

where 

9r r 2(r • v) 9r 

I - T i 2[l x (£ x l)] ®13) 

r 3v H 

and 

£2 = \ x z) x £] + Ax < e14 ) 

r 9v h l j 9v 

In these equations, Vpj is the down-range speed, and the required partial derivatives 
are computed with the following relations. 


and 



2ar 

u 


2 

r 


av H 2 /v 2 r 

6 V 



(El 5) 


8r 


3v 


2 2|j.e 


H 


(E16) 



ar 
7T 

M- 



(E17) 


where \x is the gravitational constant times the mass of the central body, and a and 
e are the semimajor axis and the eccentricity of the conic, respectively. If the 
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procedure outlined for computing flight -path -angle errors is used, equations similar 
to equations (E10) and (Ell) can be written as 


(o \ 2 = £ (t) T E (t)y_ (t) (E 18) 

\ 7TU/ 


for the mean-squared uncertainty in periapsis radius, and as 


(°r ^ = y(t) T X(t) Z (t) (E 19) 


for the mean-squared dispersion in periapsis radius. 

The second method which can be used to compute the periapsis -radius errors re- 
quires the computation of the state -transition matrix between an arbitrary time t 
along the conic and the nominal time of periapsis passage T . For example, consider 

the computation of the radius -of -periapsis dispersion. A dispersion matrix at is 

computed from 


X 




(E20) 


and, if represents the transformation from a planet- centered inertial coordinate 

system to alocallylevel coordinate system (described in the text of this paper), then 


: '{\)- A C T ,) x N A {\) 


(E21) 


from which the rms radius -of -periapsis dispersion is obtained from 




If E(t) replaces X(t) in equation (E20), the rms radius -of -periapsis uncertainty can 
be computed in a similar fashion. 
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